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Abstract 

We compute the off-forward diagonal (non-skewed) non-singlet generalized parton 
distribution of the pion in two distinct chiral quark models: the Nambu- Jona-Lasinio 
model with the Pauli-Villars regulator and the Spectral Quark Model. The analysis 
is carried out in the impact-parameter space. Leading-order perturbative QCD evo- 
lution is carried out via the inverse Mellin transform in the index space. The model 
predictions agree very reasonably with the recent results from transverse-lattice cal- 
culations, reproducing qualitatively both the Bjorken-x and the impact-parameter 
dependence of the data. 
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1 Introduction 



Recently, transverse-lattice calculations have provided first data [1] on the 
impact-parameter dependent diagonal (non-skewed) non-singlet generalized par- 
ton distributions of the pion. Generalized parton distribution (GPD) have been 
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a subject of intense studies in recent years [2,3,4,5,6,7] (for a review see, e.g., 
Ref. [8,9] and references therein) providing a unified framework for numerous 
high-energy phenomena. The impact-parameter-space formulation has been 
pursued in Refs. [10,11,12,13,14]. Actually, this is the natural framework for 
the transverse lattice QCD formulation [15,16,17,1]. In addition, the diagonal 
(non-skewed) distributions incorporate radiative corrections according to the 
standard DGLAP evolution equations for not-too-small values of the impact 
parameter b [18,19]. The results of Ref. [1] may also provide some guidance 
on the yet unknown low-6 evolution of the GPD's. 

In this paper we obtain theoretical predictions for the GPD from two differ- 
ent chiral quark models, i.e., models incorporating the dynamical chiral sym- 
metry breaking: the recently-proposed Spectral Quark Model [20,21] and the 
Nambu-Jona-Lasinio model with the Pauli-Villars regulator [22,23,24,25,26]. 
For these models it has already been shown that the 6-integrated (forward) 
parton distribution functions agree remarkably well with the phenomenolog- 
ical parameterization at Q 2 = 4 GeV 2 [27]. Our very simple predictions for 
the GPD, pertaining to a low scale of about 320 MeV, are then evolved with 
the help of the standard DGLAP equations to the scales corresponding to 
the transverse-lattice calculations [15,16,17,1]. After the evolution the results 
of Sec. 6 are in a good qualitative agreement with the data, showing similar 
Bjorken-x dependence in the corresponding impact-parameter bins. 



2 Definitions 

The off-forward (Aj_ ^ 0) diagonal (£ = 0) generalized parton distribution of 
the pion is defined by [12] 1 

H(x, f = 0, -A 2 ± ) = j d 2 b j ^l e ^-+Ax-b) 

X (n^pi)\q(0,- Z Y ,hh+ q (0, Z Y ,b)\n+(p)), (1) 

where x is the Bjorken x, and = p' — p lies in the transverse plane. This 
function has the interesting properties, 

i 

J dxH(x, 0, -Ai) = F(-A 2 ± ), H(x, 0, -Ai = 0) = q(x), (2) 
o 

1 We drop the quark flavour index since, e.g., for a positively charged pion, tt + , one 
has H u (x, 0, t) = H 3 (l - x, 0, t). 
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Fig. 1. The diagram for the evaluation of the generalized parton distribution of the 
pion in chiral quark models. 

relating it to the pion electromagnetic form factor, F(t), and to the pion 
forward parton distribution, q(x). One can introduce the impact-parameter 
representation [12], 



?(b,*) = J^e-^H(x,0,-Al)= ]^±J (bA ± )H(x,0,-Al).(3) 

where the cylindrical symmetry has been used. The second of Eq. (2) corre- 
sponds to J d 2 bq(x, b) = q(x). 



3 Evaluation in chiral quark models 



In chiral quark models the evaluation of H at the leading- N c (one-loop) level 
amounts to the calculation of the diagram of Fig. 1, where the solid line denotes 
the propagator of the quark of mass uj. Formally 2 , this yields 



H(x,0,-Al;cu) 



d A k 



rTr 



ft J (27T) 

x 5 k + - (1 - x)p ] 
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75" 



^ — ]}' — uj 1/k — uj ft — j> — uj 



(4) 



with f n = 93 MeV denoting the pion decay constant and p 2 = p' 2 = m 2 . The 
light-cone coordinates are defined as 



k + = k° + k 3 , k~ = k° - k 3 , k ± = (k 1 , k 2 ), dk°dk 3 = -dk + dk~. (5) 

2 The gauge invariant regularizations, allowing to shift the momentum in the inte- 
gral, will be specified later. 
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The calculation becomes simplest in the Breit frame, A + = 0. The Cauchy 
theorem can be applied for the k~ integration [28], yielding, after integration 
and in the subsequent chiral limit of — > 0, the result 



H(x,0,-A 



N c u 2 rd 2 K 



u = 



1 



1 + 



K ± -A ± (l-x) 



tt/2 7 (27r) 2 (K ± + (l-:r)A ± ) 2 + cu 2 



(6) 



where the relative perpendicular momentum is Kj_ = (1 — x)p± — xk±. 

To proceed further, we need to specify the regularization. First, we consider 
the recently proposed Spectral Quark Model [20,21]. The approach is successful 
in describing both the low- and high-energy phenomenology of the pion, and it 
complies to the chiral symmetry, including the anomalies. The model amounts 
to supplying the quark loop with an integral over uj weighted by a quark 
spectral density p(u>), 



H(x, 0, -A 2 ± ) = j ' duj p(uj)H(x, 0, — A^_; to), (7) 
c 

where C is a suitably chosen integration contour in the complex u space [21]. 
Next, we apply the simple techniques described in detail in Ref. [21], use the 
Feynman trick for the two denominators in Eq. (6), and integrate over K^. 
The result is 



N r \ (1 — r"l 2 A 2 

H(x, 0, -Ai) = 1 + — ^- / ou 2 p(u)du / da— ^ ' ± x , A9x .(8) 

V ' ' ±J Sn 2 f 2 J yK J J u 2 + a(l-a)(l-x) 2 Ai) K ' 

Note the correct normalization condition, F(0) = 1. Moreover, the pion elec- 
tromagnetic ms radius is (r 2 ) = — 6dF (t) / dt\ t=0 = N c /(4ir 2 f 2 ). 

In the Meson Dominance variant [21] of the Spectral Quark Model the relevant 
part of the spectral function has the form 



( , 1 37r 2 mffi 1 1 
PV[UJ) 2m AN C u{m 2 p /A-u 2 f/ 2 ' [) 

where m p = 770 MeV is the mass of the p meson 3 . The function pviyj) has 
a single pole at the origin and branch cuts starting at ±m p /2. The contour 
C encircles the branch cuts, i.e., starts at — oo + iO, goes around the branch 
point at —m p /2, and returns to — oo — iO, with the other section obtained by 



3 In this case the relation m 2 = 24ir 2 f%/N c holds [21]. 
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a reflexion with respect to the origin [21]. In the Meson Dominance model we 
get from (6) and (9) the explicit result of an appealing simplicity, namely 



m 2 (m 2 - (1 - x) 2 A\) 

H( X , 0, -A 2 ± = P \ P ^ \ 2 A 2 Y2 10 

(m 2 + (1 - x) 2 A 2 ± ) 2 



We check that H(x, 0, 0) = 1 [21] and /„* dxH(x, 0, t) = m 2 /(m 2 + 1), Eq. (2), 
which is the built-in vector-meson dominance principle. We pass to the impact- 
parameter space by the Fourier-Bessel transformation (3) and get 



q(b,x)- 



2n(l- X y 



1 — x I 1 — x \ 1 — X 



(11) 



In the Nambu-Jona-Lasinio model with the Pauli-Villars regularization one 
can proceed along similar lines as above to get 



H(x,0,-A 2 ± ) = l-^Y.^Jd* 



N r M 2 ^ } , (I- x) 2 A 2 



8n 2 f 2 Y l { M 2 + A 2 + a{l-a){l-x) 2 A 2 ± ) 

N c M 2 (l-x)\A L \ ^ fs t + (l-x)\A ± \\ 
1 + 7T272T E c * lo § 



An 2 f 2 s t ^ *\s t -(l-x)\A_ 



Si = y[{\ -x) 2 A\ + AM 2 + AA 2 , (12) 

where M is the constituent quark mass, A, are the PV regulators, and q are 
suitable constants. For the twice-subtracted case, explored below, one has, for 
any regulated function F, the operational definition [25] 



£ c t F(A 2 ) = F(0) - F(A 2 ) + A 2 dF(A 2 )/dA 2 . (13) 

i 

In what follows we use M = 280 MeV and A = 871 MeV, which yields 
U = 93.3 MeV [25]. 

It is interesting to notice that, quite generally, the chiral quark model re- 
sults displayed above depend on the momentum Aj_ and x only through the 
combination (1 — x) 2 A\. Consequently, in the b space they depend on the 
combination b 2 /(l — x) 2 . Due to this property we have 

jd 2 bb 2n q(b,x) _ /u2n ^ 2n/u2n 



J d 2 bq(b, x) 



{b 2n )(x) = (l-x) 2n (b 2n }(0). (14) 
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This means, that all the moments except for n = vanish as x — > 1, or, in 
other words, the function becomes an infinitely-narrow S function in this limit. 
This general prediction of chiral quark models is clearly seen in the lattice data 
of Ref. [1], cf. Fig. 2(b). 



4 Smearing over b 

Our aim is to compare our results, after a suitable QCD evolution, to the 
transverse-lattice data of Ref. [1]. These data give the non-singlet diagonal 
parton distribution of the pion at discrete values of the impact parameter b, 
corresponding to a square lattice with spacing of 6 ~ 2/3 fm. It is certainly 
not obvious how to compare discrete data to a continuum model. Clearly, we 
cannot achieve the continuum limit on transverse lattices, on the other hand 
we do not intend, in a simple study as presented here, to put chiral models 
on the lattice. A simple and reasonable comparison [29] is expected when the 
model predictions are smeared over square plaquettes, the same ones as in the 
discrete lattice. The plaquettes are labeled which means that they are 
centered at coordinates (ibo,jbo), and have the edge of length bo = 2/3 fm [1]. 
The smeared GPD is defined as 



(i+l/2)6o 0'+l/2)6 

V(x,[i,j]) = J dh J db 2 V(x, ^bl + b 2 2 ). (15) 

(i-l/2)6o (j-l/2)6o 

Figure 2 shows the results of this smearing. In addition, the degeneracy factor 
of the number of plaquettes equidistant from the origin is included, i.e., the 
[1, 0], [1, 1], and [2, 0] plaquettes are multiplied by a factor of four, while [2, 1] 
would be multiplied by eight. 

We note that the smearing has a large effect for the [0, 0] plaquette. This is 
because in the limit of x — > 1 the function V(x, b) becomes a distribution in 
b, which can be seen immediately from the explicit form of Eq. (11). Thus, 
the results for [0, 0] are sensitive to the size of bo- For lower values of bo the 
function becomes very sharply peaked at x — 1. 

Figure 2(b) shows the data from the transverse-lattice calculations shown by 
Dalley in Ref. [1] . These data correspond to the scale Q ~ 500 MeV, as inferred 
in Ref. [17] from the analysis of the pion light-cone wave function. Since the 
scale pertaining to our calculation is much lower, we need to evolve our results 
upward before comparing to the data of Fig. 2(b). 
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chiral models, Q=Q 



(a) 




Fig. 2. Valence impact-parameter dependent diagonal GPD of the pion, V(x, b), 
plotted as a function of the Bjorken x variable, (a) The results of the chiral 
quark models at the model scale of Q = Qo = 313 MeV. Solid lines: the Spec- 
tral Quark Model of Ref. [21], dashed lines: the Nambu-Jona-Lasinio model with 
two Pauli-Villars subtractions. Label all denotes the forward distribution, i.e., the 
function V(x, b) integrated over the whole 6-plane. Labels denote the function 
V(x,b) integrated over the square plaquettes centered at coordinates (ibo,jbo) of 
the edge of length bo, times the degeneracy of the plaquette (see the text for details). 
Following Ref. [1], the value of &o is taken to be 2/3 fm. (b) The results for V(x,b) 
at the scale Q ~ 500 MeV, obtained from transverse-lattice calculation of Ref. [1]. 
Labels as in (a). The model results of (a) can be compared to the data of (b) only 
after a suitable QCD evolution. 



5 QCD evolution 



The simple calculation of Sec. 3 has produced distributions corresponding to 
a low quark model scale, Qq. A priori, the value of Qo is not known. The way 
to estimate it is to run the QCD evolution upward from various scales Qo 
up to a scale Q where the data can be used. Alternatively, one may use the 
momentum fraction carried by the quarks at the scale Q and the downward 
QCD evolution in order to estimate Q [22,24,25]. We use the LO evolution 
with 
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/4tt\ 1 

a(Q) = UJ icgwvAW (16) 

where /3 = HC A /3 — 2N F /3, Ca = 3, and N F = 3 is the number of active 
flavors. We take A QCD = 226 MeV, which for Q = 2 GeV yields a = 0.32 [30]. 
Then one proceeds as follows: The valence contribution to the energy momen- 
tum tensor evolves as the first rr-moment of the valence quark distribution, 



y.W) ( mo) Y ?s,i2M lm 

VLWo) \a(Q )J • ( ' 

where •yf s /(2p ) = 32/81. The value of V\{Q) has been extracted from the 
analysis of high-energy experiments. In Ref. [27] it was found that at Q = 2 GeV 
the valence quarks carry 47% of the total momentum of the pion, e.g., for n + 



V x = (x (u^ - u n + d n - d n )) = 0.47 ± 0.02 at Q = 2 GeV. (18) 
The downward LO DGLAP evolution yields at the scale Q 



V 1 (Q ) = 1, G 1 (Q ) + S 1 (Qo)=0. (19) 

with G\ and Si the gluon and sea momentum fractions, respectively. The 
scale Qo defined with this prescription is called the quark model point, since 
obviously in effective quark models all the momentum is carried by the quarks. 
At LO the scale turns out to be [22] 



Q Q = 313+^0 MeV. (20) 

This is admittedly a rather low scale, but one can still hope that the typical 
expansion parameter a(Q ) / (2n) ~ 0.34±0.04 makes the perturbation theory 
meaningful. Actually, the NLO analysis of Ref. [24] supports this assumption. 
In addition, this is the same scale used in Ref. [31] to compute the pion LC 
wave function. 4 



Following Ref. [18], we apply the DGLAP evolution to the off- forward diagonal 
distribution function with the evolution kernel that does not depend on Aj_, 
or, in the impact-parameter space, on b. Then, at LO the DGLAP evolution 
in the index space simply reads 

4 An analogous analysis applied to the data of Ref. [1] shows that the momentum 
fraction carried by the valence quarks is 72% [29], which at LO would imply the 
scale of 477 MeV, compatible with the scale of 500 MeV quoted by the authors of 
Ref. [1]. 
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V n {Q,b) = J dxx n V{x,Q,b) = I ^ 



a(Q) 



7^ S /(2/3o) I 



(Qo), 



J dxx n V(x,Q ,b), (21) 



where the anomalous dimension is 



In =-20* 



3 2 4V 1 

. + (n + l)(n + 2)~ fc ^fcJ 



(22) 



with Oi? = 4/3. With n treated as a complex number, which requires an 
analytic continuation of both V n (Qo,b) and 7^ s , Eq. (21) can be inverted 
using the inverse Mellin transform 



V(x,Q,b)= J ^x- n - l V n (Q,b)- 



(23) 



The procedure, carried out numerically, is fast and stable. Since the singularity 
structure of V n (Q, b) is the same as for the forward case, we may use the 
standard Mellin integration contour in Eq. (23). 

An interesting feature of the above evolution is the induced suppression at 
x — > 1. Thus, using known methods from the b— integrated case [32], a function 
which originally behaves as V(x, Q , b) — > C(b)(l — x) N evolves into 



4Cp . a(Q) 

V(x,Q,b) -> C(6)(l - x) N ~— Xog ^) , 



x — > 1. 



(24) 



In order to compare to the transverse-lattice data of Ref. [1], we apply the 
evolution to the smeared functions of Eq. (15). Thus, we have explicitly 



+ioo 



V(x,Q,[i,j\) 



— x - n ~ l ( a ^ 
2ni \a(Q ) 



7^ S /(2A)) 



dyy n V{y,Q 0l [i )3 }), (25) 



where the distribution at the scale Q is the prediction of either of the two 
considered chiral quark models. 

We also note that in the Spectral Quark Model 



V n (Q ,b) = 

m 2 p T(n + 1) 
7r2 n + 3 



(26) 



'b 2 m 2 



bm p G 2A 



n— 1 n 
2 ' 2 



4 _I _I I 

L i 11 2 ' 2 



G 4,0 ( h2m l 



r 2,4 



n n+l 
2 ' 2 

4,o,o,o, 
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where G denotes the Meijer G function. This form can be useful for further 
analytic considerations. 



6 Results and conclusions 

Figure 2 (a) shows the plaquette-averaged functions V(x,Q ,[i, j]) for the 
Spectral Quark Model (solid lines) and the NJL model (dashed lines). We 
note that the predictions of the two models are qualitatively the same, with 
the NJL curves pushed to somewhat lower values of x. For the lack of space, 
in this paper we display the QCD evolution of the Spectral Quark Model only. 
The case of the NJL model is qualitatively the same, with the corresponding 
curves moved to a bit lower values of x, simply reflecting the different initial 
condition of Fig. 2 (a). These results and other details will be presented in a 
longer paper. 

The results of the evolution are shown in Fig. 3 at three values of the reference 
scale Q: 400 MeV (a), 500 MeV (b), and 2 GeV (c). We note a large effect of 
the evolution on the distribution functions. The lines labeled all correspond to 
the forward case, i.e., show / d 2 bV(x, Q, b) = V(x, Q, Aj_ = 0). The originally 
flat distribution of Fig. 2(a) recovers the correct end-point behavior at x — > 1 
according to Eq. (24). As Q increases, the distribution is pushed towards lower 
values of well known for the DGLAP evolution. At Q = 2 GeV the 

result agrees very well with the SMRS parameterization of the pion structure 
function [27], as can be seen from Fig. 3(d) (here we plot for convenience 
xV(x,Q)) by comparing the dashed and solid lines. This result was already 
obtained in Refs. [22,24]. 

The results for the plaquette [0, 0] follow, at large x, the forward distributions. 
This is clear from the behavior described at the end of Sec. 3, i.e., from the 
dependence of the initial function on the variable bj (1—x). Certainly, as x — > 1, 
the integration over the [0, 0] plaquette is the same as the integration over the 
whole 6-space. At Q = 400 MeV and 500 MeV the values of V(x,Q, [0,0]) 
reach a maximum at an intermediate value of x, and develop a dip at low x. 
This is in qualitative agreement with the transverse-lattice data of Fig. 2(b). 
We note that there the dip at low x is lower than in our model calculation, yet, 
in view of the simple nature of our model and approximations (chiral limit, 
LO evolution, evolution independent of b, uncertainties in the determination 
of b and Q on the lattice) the similarity is quite satisfactory. We have checked 
that if the value the lattice-spacing parameter, bo, were lowered, an even more 
quantitative agreement would follow. 

The results for non-central plaquettes also qualitatively agree with the lattice 
measurements. In this case at x — > 1 the corresponding functions vanish very 
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Fig. 3. Results of the LO DGLAP evolution of the impact-parameter depen- 
dent diagonal non-singlet generalized parton distribution function of the pion, 
V(x, b, [i, j]), started from the initial condition at Q = Qq = 313 MeV produced 
by the Spectral Quark Model (Fig. 2(a), solid lines). Figures (a,b,c) correspond to 
Q = 400 MeV, 500 MeV, and 2 GeV, respectively. Labels as in Fig. 2. Figure (d) 
shows xV(x,b, for Q = 2 GeV, with the dashed line showing the SMRS [27] 
parameterization of the data for the forward parton distribution function. 



fast, in accordance to our model formulas. The difference with the lattice 
calculation of Fig. 2(b) is that in our case the farther plaquettes naturally 
bring less and less, and the yield from the [2, 0] plaquette is lower than for the 
[1, 1] plaquette. In Fig. 2(b) it is the other way around. 



In summary, the obtained agreement of our approach, based on non pertur- 
bative chiral quark models in conjunction with perturbative LO DGLAP evo- 
lution, with the data from the transverse lattices, is quite remarkable and 
encouraging, baring in mind the simplicity of the models and the apparently 
radically different handling of chiral symmetry in both approaches. We also 
note that the low-energy scale taken for the chiral quark models is consis- 
tent with previous analysis based both on the forward parton distribution 
amplitudes as well as the light cone wave function. Our analysis might be 
reinforced by extending our calculation to include the NLO perturbative cor- 
rections. Such a study is left for a future research. 
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